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Abstract. A statistical mechanics model of isotropic turbulence that renormalizes the effects of turbulent 
stresses into a velocity-gradient-dependent random force term is introduced. The model is well-defined 
within the context of the renormalization group e expansion, as the effective expansion parameter is 0(e). 
The Kolmogorov constant and N parameter of turbulence are of order unity, in accord with experimental 
results. Nontrivial intermittency corrections to the single-time structure functions are calculated as a 
controlled expansion in e. 

PACS. 47.27. Ak Fundamentals - 47.27.Gs Isotropic turbulence; homogeneous turbulence - 05.10.Cc 
Renormalization group methods 

Key words. Renormalization group, turbulence, inter- has focused upon the random force model that Domini- 
mittency, operator product expansion cis and Martin generalized from Forster, Nelson, and 

Stephen's model of a randomly stirred equilibrium fluid 
Pj- This model has received quite a bit of attention, and 

1 Introduction 

its transport properties have been further examined by 
Turbulence remains an outstanding problem in classical Yakhot and Orszag and Avellaneda and Majda and co- 
mechanics. Early on, several self-consistent or closure-based workers |4ll51lrj] . Various treatments, up to and including 
statistical models were introduced, with perhaps Kraich- field-theoretic e expansions, have been applied to the ran- 
nan's Direct Interaction Approximation the most well known dom force model of turbulence. As pointed out by Eyink, 
P]. Within the statistical mechanics literature, attention however, even the e expansion does not lead to a controlled 
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calculation in the random force model because the term that are quenched in time. Moreover, in practical, engineering 
representing the effects of high Reynolds numbers must be type calculations, the turbulent forces are often related to 
calculated to arbitrarily high order in perturbation theory, velocity gradients by a constitutive relation that contains 
even for an 0(e) calculation More recently, a variety of a "turbulent viscosity" parameter, the simplest of these re- 
novel renormalization group techniques has been applied lations being linear |13|. Following this line of reasoning, 
to the problem of Navier-Stokes turbulence |XI9ll()llf| . we introduce a new statistical mechanics mode] for tur- 

m, . , ii i i i-ii- c . bulcncc that rcnormalizcs the effect of turbulent stresses 

the experimentally observed scaling behavior of tur- 

,, ,. . c, ,, , ,i T , , into a velocity-gradient-dependent term. This model will 

bulent energy dissipation, often called the Kolmogorov en- J ° r 

, ,, , n i turn out to be well-defined within the e expansion. That 

ergy cascade, suggests that there may be a strong analogy 

i , , , , , , T , , , n is, the renormalization group theory takes into account all 

between critical phenomena and turbulence, tndeed, the ° 

i r , , , , , e ,i n , , physical effects of this model, consistently to order e. Due 

search tor such an analogy motivated much ot the held- 

.-, , . i n . , , , , rm t ^1 to its similarity with turbulent viscosity type models, this 

theoretic work on fluid turbulence [12]. Just as there are J J ' 

., , n -ii i • ■ approach may lead to a closer connection with practical 

wild fluctuations m particle density near a critical point, rr J ^ 

, n • ,i a -j i , turbulence calculations, 

so to are there large fluctuations m the fluid velocity at 

high Reynolds number turbulence. This similarity sug- 
gests that the effects of turbulence may be modeled by 
a random, velocity-dependent force, just as the effects of 
critical fluctuations can be modeled by a random, density- 
dependent force in the standard 4> A model. While there 

are likely random forces that arc independent of the fluid TTr . , , , ,. , , , n 

We introduce our velocity-gradicnt-dcpcndcnt random 

velocity in the context of turbulence, there are also very r , , . „ „ m, n n , . ,, r n 

force modef in bee. I. the probfem is cast m the framework 

likely random forces that are velocity dependent. Thus, c , ■ i j , n u ,i on t 

J 'of time-dependent held theory m bee. 6. Renormalization 

the turbulent force should really depend on the gradient a . . , , , , , . 

J 1 ° group flow equations are also calculated in this section, 

of the velocity, as it is not simply large velocities that lead , . , n . r , , n . 

and three appendices provide details of the renormahza- 

to turbulence, but rather regions of high gradient, such as . , ,. m u , • ri i a 

' b tion group calculations, the behavior of these flow equa- 

walls or boundaries, that lead to turbulent behavior. Such . . , , . , , ,. . 

' tions m two and three dimensions, and the predictions 

boundary roughness is one mechanism for breaking the r T ^ , , n ., , . 

J ° ° for the Kolmogorov energy cascade, are described m bee. 

symmetry from laminar to turbulence fluid flow. . t> t , . . , . , ., , ,. , , . 

J J 4. JNontnvial mtermittency corrections to the single- time 

Indeed, random boundary roughness along the walls structure functions are calculated by an operator product 

generates velocity-gradient-dependent forces in the bulk expansion in Sec. 5. We conclude in Sec. 6. 
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2 Velocity-Gradient-Dependant Random Initially, we treat this as a mathematical problem, tak- 

Force Model ^ ^° ^ e ar ^ rar y- Later, we determine y by requir- 

ing that the transport properties of turbulence are repro- 

Our goal is to write a form of the Navier-Stokes equation 

duced. This scaling form of the correlation function, Eq. 

that contains a random piece, the random piece represent- 

l|3|l. applies only in the incrtial, Kolmogorov regime, for 

ing the statistical effects of the turbulence. The Navier- 

wavevectors below an upper cutoff related to the inverse 

Stokes equation with a random force is 

of the dissipation length scale and above a lower cutoff re- 

dtVi + Y\ ik J2j VjdjVk = uV 2 Vi + ft , (1) lated to the inverse of the so-called integral length scale. 

It is this Kolmogorov scaling regime that is of interest in 
where v = p/ p is the kinematic viscosity, and /; = IIi k (F k — 



the present work. Given the form of Eq. [21 the parameters 
\f~D a and can be viewed as modeling gradients of 

the turbulent viscosity. 



dkP)/p is the total body force on the fluid. The pres- 
ence of the projection operator IIik(k) = 8i k — kikk/k 2 
in these formulas ensures that the incompressibility con- 
dition V • v = is maintained [3]. The Fourier transform 

is defined by /(k) = J tfx/(x) exp(ik x). 3 Renormalization Group Calculations 

We choose the random force to depend on the gradient 

c , , , . We write the Navier-Stokes equation in field-theoretic form 

of the velocity: H 



/i(x,t) = -j ijk (x)djVk(-x.,t) . (2) 



so that the renormalization group can be applied sys- 
tematically within the e expansion |TJ. Within the field- 
We assume that -y ijk is symmetric under exchange of j and theoretic formalism, any observable can be calculated. The 
fc, so that the force looks like a turbulent stress. We aver- average velocity, for example, is given by v l (x, t) = (6, (x, t)) , 
age over the statistics of this force using a field-theoretic w here the average over the b field is taken with respect to 
representation of the Navier-Stokes equation. We choose the weight exp(-S). Using Eqs.UHl we arrive at the fol- 
the correlation function to be lowing action: 

(7 i , fe (k 1 )7^(k 2 )) = (2 7 r)^(k 1 + k 2 )|k ir ^X. (3) g= j [ d a i{ - k , m + uk * + 5(t)MKt) 



where 



f f dt (2 7 r) d <5(k 1 + k 2 +k 3 ) 



D l ™ k n = D a {8 a 8 jm 5 kn + 5 in 6 jm S k i x fc^S* (k u t)bi(k 2 ,t)bf(k 3 ,t) 



+ Sii5j n 5km + SimSjiSkn + Si m 6j n Skl + Si n Sjl5km) 



If [ dhdt 2 (27r) d 8(k 1 + k 2 + k 3 + k 4 

1 Jkik 2 k 3 k 4 J 



-Dp^ijSlmSkn + SikSjmSln + StkSjnS^ + 8ij8 km 8i n ) . X k 2j k4 m D ^ k |ki + k 2 1 v 

(4) X ^(k 1 ,t 1 )6^(k 2 ,t 1 )6 / " L (k 3 ,t2)^(k4,t2) 
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(5) 

The notation J k stands for J dk/(2-7r) d , the integrals over 
time are from t — to some large time t = tf, and the 
summation convention is implied. This action is written 
in terms of the divergence- free part of the field, bf- (k) = 
J2k ^ifc(k)&fc(k). We have used the Feynman gauge, adding 
in a curl-free component in the quadratic terms to make 

Fig. 1. (a) Diagram representing the propagator. The arrow 

later calculations easier. Initially A = 1. We have used the 

points in the direction of increasing time, and double lines rep- 
replica trick |15| to incorporate the statistical disorder, 

resent the bar fields, (b) Convection vertex A. (c) Disorder 

but have suppressed these details since they do not enter 

vertex D a and Dp. 

in a one-loop calculation. 

We now apply renormalization group theory to this ac- <Ai±. The b field is scaled as 6'(k',i') — a z 1 d+a b(k,t). 

tion. It is important to note that the fields must all have To maintain dimensional consistency so that the b field 

zero average value before the renormalization group is ap- scales as a velocity, one must set a = [17] . To keep the 

plied Q21[n|, otherwise it would not be correct to truncate tim e derivative in S constant 18 , the b field is scaled as 

perturbation theory at any finite order. If there were an b ( k '> O = a 1 -^ Q 6(k, t). In the loop calculation, we make 

average velocity, the action in © would be different, con- use of the relation for the reference system averages 

taining a term of the form — ik ■ (v) in the propagator. , ,, , , . 

(6( (k 1 ,t 1 )6 J ± (k 2 ,t 2 )) = (27r) d 5(k 1 +k 2 ) 

The vertices in the theory are shown in Fig. ^ 

x77y(k 2 )Go(k 2 ,i 2 -ii) 

From power counting, the upper critical dimension for 

(6) 

this theory is d c — 2 + y. Note that this upper critical 

dimension is exactly defined once the model is specified where Go(k, t) — exp[— vk 2 t]0(t), and 0(t) = 1 if t > 
|14| . The deviation of the physical dimension from the and otherwise. Note that elimination of modes at one 
upper critical dimension is parameterized by e = d c — d. end of the spectrum by perturbation theory is the stan- 
We use the momentum shell procedure, where fields on dard procedure in renormalization group theory |14|. Use 
a shell of differential width In a = dl are integrated out, of Eq. © does not imply that the system is somehow 
A/a < k < A. Note that the combination dl invariably Gaussian, as the parameters within the renormalized t he- 
means a differential on I; in all other cases, the factor d de- ory are flowing. The critical properties of the Ising model 
notes the physical dimension. As usual, we rescale time by at the non-Gaussian Wilson-Fisher fixed point, for exam- 
the dynamical exponent if = a~ z t and distance by k' ± = pie, are analyzed in exactly this way |14j . We make use of 
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lar diagrams of the form Fig. [2b, detailed calculation of 
which is described in Appendix B. Finally, there are 100 
diagrams of the form Fig. and a detailed discussion of 
them is given in Appendix C. These last diagrams require 
some care in their calculation, as they contribute to the 
complex tensor structure of D^J 1 in Eq. JSJ. These last 
diagrams make contributions in the form 

k 2j k 4v b i (ki,ti)6jt (k 2 ,ii)6 r {k 3 ,t 2 )b w (k 4 ,t 2 ) 



v 7i fstu T~\lmn t~\uvw 
x m \mn U ijk u rst 



(8) 



Fig. 2. One-loop diagrams: (a) self-energy diagrams contribut- 
ing to v, (b) vertex diagrams contributing to A, and (c) vertex 
diagrams contributing to D a and Dp. 

the rotational averages: {kik u )n = 5i u k 2 /d, (kik m k s k u ) q = 
{5i m 8 us +8i u 5 ms +8i s 5 mu )k A /[d(d+2)], and (h 
M^ n k 6 /[d(d+2)(d+4)], where the function M$? n is equal 
to all possible couplings of pairs of the arguments: 

Mimn = 3lt8ms8nu + 5lt&mn5 S u + Sl t S mu 5 ns 
+ 5l m 8t s 5 nu + Si m 5 n tS su + SlmStuSns 
+ &ls&mt&nu + dlsSmnStu + Sl s 5 mu S n t 
+ $ln$mt&su + 5l n 8 ms 5tu + Sl n 5 mu S s t 

(7) 

The one loop contributions arc shown in Fig. El There 
are 13 diagrams of the form Fig. detailed calculation 
of which is described in Appendix A. There are 13 simi- 



There are 15 terms of the form Fig. when the sums 
over I, m, n and s, t, u are taken. Each of these terms cor- 
responds to one of the 15 terms in Eq. Q, and each gives 
a contribution that is exactly of the form in Eq. JSJ. The 
theory is, therefore, self-consistent in that no terms are 
generated at this order that are not of the original form, 
and the symmetry of the original theory is maintained. 

The flow equations that result from the one-loop cal- 
culation arc 



<iln v 
dl 



dlnA 
dl 



z-2 + [{d 2 + 2d - 2)D a + {d- 2) Dp] 



yK d 



v 2 d{d + 2) 



A~ 



-a - [(d 2 + 2d- 2)D a + (d- 2)Dp] 



-A- 



dhxD D 
dl 



v 2 d{d + 2)' 
2z - 4 + e - [2D a + 4D P ] 

2K d (d 2 + 2d-2) 



dlnD 



dl 



v 2 d(d + 2)(d + 4) 
2z - 4 + e - [2(d 2 + 2d- 2)E^/D 

+ (d 3 + 6d 2 - 8)D a + {d 2 + 2d- 8) Dp] 
2K d 



v 2 d(d + 2)(d + 4) 



A- 



(9) 
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The constant K d = S , d /(2 7 r) d , and S d = 2n d ^ / T(d/2). mimic the effects of, say, wall roughness, break statisti- 
In solving these equations, we set a — for dimensional cal Galilean invariance, and this allows X(l) to flow, in 
consistency. In the standard model of turbulence |2], terms contrast to the conventional model with random forces 
higher order in X(l) must be kept in the flow equations. In delta-correlated in time. Indeed, our explicit calculation 
the present model, A(/) flows to zero rapidly, and higher shows that X(l) decays exponentially to zero, 
order do not contribute at this level. We now turn to a calculation of the energy spectrum, 

defined by PU 

4 Results and Discussion E(k) = ^ - 1 ^ Kdk^C^k) , (If) 

In two dimensions, both D a and Dp are relevant. We find where the velocity- velocity correlation function is given in 
a fixed point of D* = D* = 2ev 2 A e /[3K 2 (l + y)}. The thc field-theoretic language as 



dynamical exponent is given by z = 2 — ye/ [2(1 + y)]. The 
Reynolds number term scales as X(l) — Xoe~ el ^ 2 ( 1+v ^ . 

In greater than two dimensions, Dp reaches a fixed 
point, but D a flows to zero. We find 



(27r)^(k 1 +k 2 )iT ij (ki)C ij -(ki) = (bf(k u t)bf(k 2 ,t)) . 

(12) 

Under the scaling of time and space that occurs within the 
renormalization group calculation, this correlation func- 
d(d + 2)v 2 A f - ti on sca i es as 



Dp 2(1 + y)(d - 2)K d 



V±_ d = > (bi'(k' llt ')bj'(k' 2 ,t')) 

2 ( 1 + y ) lJ a d (2w) d S(k' 1 +k^n ij (k[) 



= 2 — C^k) = : 

A(l) - Xoe-WW+v)] . (io) = ad+ 2-2 Z -2 ad , {W) (13) 



Interestingly, the dynamical exponent is the same in two Making the assumption that Cy(k) - (const) k~ s , we find 

and greater dimensions, as is the decay of X(l). f rom Ea.HSIthat 5 = d+2-2z-2a. The energy spectrum, 

A particularly beautiful feature of this theory is that therefore, scales as 
X(l) decays exponentially to zero. This property is what 

E(k) ~ (const)fc 2z - 3+2Q . (14) 

makes the theory well defined within the epsilon expan- 
sion. If A(Z) had stayed at unity, the vertex in Fig. \l]p could For any isotropic statistical theory of turbulence, then, 



be inserted arbitrarily many times in the loop expansion, the dimensional consistency condition of a = enforces 

and terminating the expansion at any finite order would a relation between the dynamical exponent and the expo- 

not be justified by any small parameter. The present cal- nent of the energy cascade. In particular, the Richardson 

culation, on the other hand, is a controlled expansion in separation law implies z = 2/3, and this result is equiva- 

e and \(l). Note that the quenched random forces, which lent to enforcing the Kolmogorov energy cascade: E(k) ~ 
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(const) fc~ 5 / 3 . The relation z = 2/3 implies y = 3.44152 in Using Eqs. (JSJl and JEjJ, we find e = (27/8)zy 3 C£/l 4 . 
two dimensions and y = 4.28849 in three dimensions. Finally, to complete the matching we assume that the 

To calculate the Kolmogorov constant Ck we introduce wavevector cutoff is one-half the Kolmogorov dissipation 
a source of randomness into the model: number, 2A = k d = (e/i/ 3 ) 1 / 4 [201 - Putting these relations 

together in three dimensions, we find 



^=-§E/ dtJj(k)bt(~-k,t)bUk,t) • (15) 



£ = 0.0380Z? 

In the range for which scaling occurs, we set c(k) <~ fc 3z ~ 2 ~ d . 

C K = 1.68 

For later convenience, we also require that \\m. r ^Q c(r) = 

■, mt , • , ■ , ■ N = 0.397 . (19) 

1. the randomness expressed m this source term drives 

the model away from the trivial solution v,-(x, t) = 0. Note AT , , n , , , , , . , 

v ' ; Note that field theory cannot calculate non- universal pa- 

that the randomness parameter scales as D(l) = D. This , , •,, , , , 

^ v ' rameters such as these with precision, as these results de- 

randomness parameter does not contribute to v since X(l) , ,i • • ™ n-n» j i ^ ^ 

v ' pend on the assumption m _bq. (|18|) and the relation be- 

flows to zero, and nothing contributes to D(l) at one loop. A A j^-ij ±. i_- i i a- 

' & v ' 1 tween A and kd. A more detailed matching calculation 

Since D(l) contributes to physical properties at higher r , , n . , , , , £ , , , ^ m 

w r^ & of these values using the model of turbulence m Fqs. (|3J) 

loops only through X(l). D n (l), and DM), all of which are , rm , n , , • , , ,. , , , r 

^ J & v " v y ' pw ' and (|15|) to refine the matching calculation would be of 

small, the effects of D are controlled within the e expan- 



mtcrcst. 



sion. Using the matching Eg. 1131 the correlation function is 
given by C tj (k) = [DA*- 2 /(2v)]k~ s . For fully-developed 
isotropic turbulence z — 2/3, and using the notation 0] 



5 Intermittency 



E(k) ~ CK£ 2 ^ 3 fc 5 ^ 3 (16) ^ e ^ eTe address the issue of intermittency in our model. 

That is, we seek to determine the scaling of the single-time 

we find C* K £ 2/3 = (d- l)K d DA-^ 3 I Uv), Similarly, the 

structure function 

wavevector-dependent viscosity considered in the fluid me- 

chanics literature [U is given by i/(k) = vA^k^' 3 . Using S ^(r) = ([Ib^x + r) - b ± (x)| 2 ]") ~ (const)/ 2 " . (20) 



the notation 0] 



From simple dimensional analysis, we find C,2n — nz. From 
v(k) ~ A^e 1 / 3 ^ -4 / 3 (17) explicit calculation for our model, we find an exponent 

that differs from this value. This difference is referred to in 
the fluid mechanics literature as intermittency. That C^n 7^ 
nz is an expression of the non-Gaussian nature of the fixed 
J o 2vk E(k) . (18) point id en tih e( i m our model and of the divergence of the 



we find Ne 1 ' 3 = vA A l 3 . The energy dissipation rate is 
given by 
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single-time structure functions in the limit of an infinite L — ► oo, and so it does not play a role. In our case, on the 
integral length scale. other hand, the scaling function gives the corrections to 

We use the same arguments about scaling of space and intermittency. 
time as in Eq. El t° express the original correlation func- To determine the function /„ we use the operator prod- 
tion in terms of the renormalized correlation function: uct expansion 14,21,22 . A similar strategy has proven 

successful in the study of turbulent transport of passive 

S2n(r)~e n ^- 2 ^ l S 2n {r{l)-l) . (21) 

scalars |23| . The operator product expansion states that 

Here r = exp(Z)r(/). This equation is applied until r(l) is of (F (x 1 ) J F(x 2 )) ~ ^c«(| Xl - x 2 |)<F«(x)) . (24) 

the order of the dissipation length scale r(l*) = 2ir/A = h. a 

Here the x = (xi + X2)/2, and is the set of all renor- 

At this length scale, we then match the correlation func- 

malized operators that are generated by the renormaliza- 

tion to a perturbation theory result: 

tion group flow of F . Equation 1241 is nothing more than 

lr 

a Taylor series expansion, where both the bare terms in 



S 2n (r(l*);l*) cx 



2^* 

, 2 \ n 

_ I „(3*— 2)ni* n (u\ the Taylor series and those terms that are generated by 

2u) { ' 

u2\ n the renormalization flow are included. In our particular 

*_\ e (3*- 2) „r f (22) 

' case, instead of a pair of operators, we have S 2n , which is 

where we have used the scaling of the small-r behavior of a product of 2n factors on the left hand side of this equa- 
the c(r) function in Eq.U^and have used lim h _ c(h) = 1. tion. The important point about this expansion is that the 
Combining Eqs. E2and[221 we find functions c a (r) are finite and exhibit no dependence on 

the system size L. Any possible system size dependence of 

S 2n (r) ~ (const)r" 2 /r l (r/ J L) . (23) 

this expansion, therefore, is contained within (F^(x.)). By 

We have here introduced the fact that traditional renor- comparison to Eq. E3 we see that the scaling function /„ 

malization group arguments can determine asymptotic be- jg thus determined by the behavior of these renormalized 

havior only up to a scaling function of the ratio r/ L, where operators. We first determine the scaling of the renormal- 

L is the macroscopic size of the system. This is because ized operators: 

r and L are scaled by the same factor in the renormal- „ . ,» 

(F«) = e A « l (F*(l*)) . (25) 

ization group analysis, and so no dependence on the ra- 

We follow the renormalization flows until L(l*) — exp(— l*)L 

tio r/L is detectable. For turbulence, h is the dissipation 

r, a criterion that automatically ensures the functional 

length scale, and L is the integral length scale. In many 

form specified in Eq. We thus conclude that 

applications of renormalization group theory to condensed 

matter systems, the scaling function f n (r/L) ~ (const) as (Fa) = (y . (26) 



Park and Deem: A Statistical Mechanics Model of Turbulence 9 

In fact, we determine the scaling of these operator averages We first determine the scaling function j\. We limit 
by using the generating functional SS = —A a J dt J d d xF^. consideration to the case d > 2, where -D* = 0. The corre- 
We will find lation function S 2n is given in a Taylor series as S^n (r) ~ 



A a (l) ~ e {z+d - 2nz+ ^ )l A a {Q) . (27) 

In terms of A a and the partition function Z , the operator 
average is given by 



[r 2 d x bid x bi] n . This will generate the symmetrized opera- 
tor F — [djbidjbi + dibjdjbi] n , which we consider. 

For the case n = 1, we consider the generating func- 
tional 



(28) 

a q (o)=o SSu = - dt 



«(0) 

This equation makes clear why Ea. l25l has the form that it 
does and identifies A a = —2nz + 7 a . The value of 7 Q will 
be determined by a nontrivial fixed point of the renormal- We nnc ^ 



A( 1 )k 1 -k 2 6 ± (k 1 ,<)-S ± (k 2 ,t) 

(32) 



kik 

A^k 1 -b ± (k 2 ,t)k 2 -b ± (k 1 ,t) 



ization group flow equations for A a . Combining Eas. l2l)l dA^ \ aU) 2(d — l)g ,^ 

— = (z + d-2z)A + - -A 

and [23 we find 



\ 2nz— , y Ci . 

j) (F?(l*)) ■ (29) 



(33) 

where g = DpKj/^A 6 ). We, thus, identify 71 = 2(d — 



l)g*/d. Using Eq. 1101 we find in three dimensions z = 
The function (F*(l*)) is determined by matching exactly 2 - e/4 + 0(e 2 ) and g* = 15e/4 + 0(e 2 ) and conclude from 
as in Eq.EU We, thus, find that Eqs.ED and EU that 



{Fa) \l) \l) S 2 (r) ~ (const)r* (£) 

f\n(2-z)-7, 



f \ 2—3 — fi 



To make use of this result of the operator product expan- 
sion, it remains only to calculate the value of 7 Q . Once we 
have this value, we find the scaling function to be 

/ r \ »(2-z)-7<» 

hn{r/L) ~ y— j as L to 00 . (31) 



= (const)r 2 - £ / 4 (£j . (34) 

We now determine the scaling function f 2 . We start 
with the symmetrized generating functional 

5S^=-fdtf b ± (k 1 ,t)-b ± (k 2 ,t) 



'kik2k3k4 

xb ± (k 3 ,t) ■ 6 ± (k 4 ,i) 



4 2) ki •k 2 k 3 -k 4 



+A {2 \k 1 • k 3 k 2 • k 4 + ki ■ k 4 k 2 • k 3 ) 



• (35) 



In fact, the operator S 2n will generate several new opera- 
tors F a in the expansion of Eq. 1241 The appropriate value 
of 7„ to use in Eq. EH is the largest one. These generated This term generates two additional generating functional: 

operators may mix upon the renormalization, in which „ „ 

5S {1) =- dt b ± (k 1 ,t)-b ± (k 2 ,t) 

of the flow equation matrix. 



case the appropriate value of j a is the largest eigenvalue •* " , kik 2 k 3 k4 

X A^ ) k 1 -k 2 k i -b ± (k 3) t)k 3 -b ± (k 4 ,t) 



10 
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and 

6SW> 



dt 



A^kx k 2i k 3i k 



kik 2 k 3 k4 



We have, therefore, derived the intermittency correc- 
tions to dimensional analysis for the present model. The 
corrections are calculated in a controlled fashion and are 
proportional to e. The coefficients of the correction are not 
small, and for finite e, higher order terms in the expan- 
sion are required for an accurate estimation of the effects 
(36) °f intermittency. 

One might wonder whether there are any corrections 
to the Kolmogorov energy cascade, Eq.^J that arise from 
the operator product expansion. More generally, are there 



(kr • k 3 k 4 • b ± (k 3 , i)k 2 • 6 x (k 4 , t) 
+k 1 -k 4 k 2 -^(k 3 ,i)k 3 -6 ± (k 4 ,i) 
+k 2 -k 3 k 4 -6 i (k 3 ,i)k 1 -6 ± (k 4 ,i) 

+k 2 -k 4 k 1 -6 i (k 3 ,;)k 3 -6 ± (k 4 ,;) 
+4 1} (k 3 ■ k 4 kx ■ ^(ka.tjka • & x (k 4l t) 
+k 3 -k 4 k 2 -^(k 3 ,t)k 1 -6 ± (k 4 ,t) 



+ A 2 ) ^ k ii k 2ik 3i k 4fe + ki J -k2fck 3 jk4 i 
+k u k 2l k 3i k 4fc -I- ki fc k 2j k 3; k 4j 
+Af ] ^ki fe k 2; k 3i k 4j + k u k 2i k 3j k 4i 
+k 1 ,k 2 fek 3i k 4 ,- +k u k 2fc k 3l k 



We will find A ( 2 1] (l) = A^\l). Although we have included 



A^ for generality, we will find that this term is not gen- 
erated, and A { ° ] (I) = 0. 



any corrections to ([6^(x, t)6^~(x, f)]™)? There are no such 
corrections. This type of operator flows under the renor- 
malization group only to operators with more derivatives, 
such as dibj-dmb^b^-b^-) 71-1 . These operators are less rel- 
evant than the original operator, and so they make no con- 
(37) tribution to the scaling at leading order. In the language 
of Eq. [21 = 1, and the scaling function f n (r/L) ~ 
(const) as L — + oo. 



A lengthy calculation shows that 



dA , g 

— = (z + d- 4z)IA + — ^ 

dl y ' d(d + 2)(d + 4) 



MA 



(38) 



6 Conclusion 

An alternative, simpler model would have been to take the 
turbulent forces to be proportional to the velocity, rather 



where the vector A = (Af' , Af ] , Ap , A$' , A^ , A^ ) , A { ° ] ) , than the velocity gradient. The simplest model, moreover, 

and the matrix M is given in Table □ We diagonalize the would take the forces to be white noise in time and uncor- 

matrix M, finding the eigenvalues (0, 0, 0, 0, d 3 + id 2 - related in each of the spatial dimensions. For this model to 

4d-16,2(rf 3 + 5rf 2 + 2d-8),4(rf 3 + 7rf 2 + 6d)). The largest be nontrivial, a mean fluid flow must be introduced [HI- 

eigenvalue is the last one for all d. In three dimensions we Interestingly, when this is done for forces that are random 

identify 7 2 = 432^*/105. Using Eqs. [221 and E2 we find i n time as well as space, the resulting theory has the same 

that fl ow equations as the traditional random force model of 

S 4 (r) ~ (const)r ~ £/ U-J . (39) turbulence HE]. Also of interest to note is that a theory 
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with velocity-gradient-dependent random forces that are 
white noise in time would have no renormalization of any 
parameter, as the diagrams of Fig. [21 would all vanish due 
to the causality of the bare propagator, Eq. JJ]J. 

The scaling of the Kolmogorov energy cascade is deter- 
mined once the value of the dynamical exponent is fixed, 
i.e. z = 2/3 for isotropic turbulence. In random force mod- 
els such as the present one, the scaling of the energy cas- 
cade simply serves to fix the correlation function of the 
random forcing. The predictive power of models such as 
these lie in their ability to provide nontrivial predictions 
of the intermittency corrections. In the present model, we 
are able to provide these corrections as a systematic ex- 
pansion in e. 

In summary, we have introduced a new statistical me- 
chanics model for isotropic turbulence. This model makes 
use of a random, velocity-gradient-dependent force. This 
model is both consistent with practical, engineering-type 
calculations and well-defined within the renormalization 
group e expansion. This model makes stronger the anal- 
ogy between turbulence and critical phenomena. 

Owing to the irrelevance of the convection terms at the 
fixed point, our results may alternatively be viewed as an 
analysis of transport in a new class of random media. 



Appendix A: One-Loop contributions to v 

We here show how the diagrams of Fig. |3l contribute to 
the propagator, Fig. Each of the terms is associated 
with one of the D l ™£ terms in Eq. JSJ. In the calculation of 
the averages on the shell, there are five terms associated 
with D a and three associated with Dp, as the last two 
terms associated with D a are identical, and the last two 
terms associated with Dp are also identical. The first term, 
associated with D a , is 
D„ 



h = 2 x 



dt\dt2 



kik 2 k 3 k 4 



x ^Trf^kx + k 2 + k 3 + k 4 )|ki + k 2 | 



X0(h - t 2 )(2Tr) d 5(k 2 + k 3 )e-^^-^ 
x jk 2 -kSi (ki,ti)6|(k 4 ,t2)^(k 2 ; 



K d 



D, 



y(d+i) 

'd(d + 2) 



dqq 



d-y-3 



A/a 



x J dtj yk% (-k,t)b±(k,t) (40) 

For the remaining contributions, we list the symmetry 
factor, terms in braces in the integrand of Eq. I|4(J|) that 
change, and the final contribution in brackets that change. 
The contributions are shown in Table |2 where the depen- 
dence of the fields upon time has been suppressed. The 
term e _I/fc ^ tl_ * 2 - ) , which has the same momentum argu- 
ment as the Tliki) term, has been suppressed. The delta 
function is also suppressed. Summing all these contribu- 
tions to we get the first flow equation of Eq. @. 



This research was supported by the National Science Foun- 
dation and by an Alfred P. Sloan Foundation Fellowship to 
M.W.D. 



Appendix B: One-Loop contributions to A 

We here show how the diagrams of Fig. |2h contribute to 
the convection term, Fig. QJd. Each of the terms is as- 
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sociated with one of the -Dfe™ terms in Eq. 0. It is delta functions have also been suppressed. Summing all 

convenient to define the convection operator My-fc(k) = these contributions to A, we get the second flow equation 

kj Tlik (k) + fcfc Tlij (k) . The convection term of Eq. (J5J then of Eq. JSJ . 
becomes 



+ % A [ [ dt (2n) d S(k 1 + k 2 + k 3 ) 

* Aik 2 k 3 J 

xM^(fc 1 )^(k 1 ,i)6|(k 2 ,t)6^(k 3 ,t) . (41) 

In the calculation of the averages on the shell, there are 
again five terms associated with D a and three associated 
with Dp. The first such term is 

ji=2x (4) » (£ ) f / ** 

x f (2^)^(k 1 +k 2 + k 3 +k 4 ) 

Jkik2k 3 k4 

x f dt 3 [ (2n) d 5(k 5 + k 6 + k 7 ) 

J A 5 k 6 k 7 

x\k 1 + k 2 \-y0(t 1 -t 3 )o(t 3 -t 2 ) 

x(2^(k 2 +k 5 )(27r) d <5(k 6 + k 3 ) 

xe -^I(t 1 -t 3 ) e -^(t3-*2) A ^ mn ( k5 ) 

x\k2 k k ik b i {k u t 1 )bf(k 47 t 2 )b^(k 7 ,t 3 ) 
xiT i /(k 2 )/7 im (k 6 ) 



2 v 2 



-D, 



(d+1) 



d— 1/— 3 



yl/a 



d(rf+2) 

x fdt f (2n) d S(k 1 +k 2 +k 3 ) 

J Jkik 2 k3 

x M ijk (k x )Si (ki , 4)6+ (k 2 , 4)6^ (k 3 , 4) 



(42) 



Appendix C: One-Loop contributions to D a , Dp 

We here show how the diagram of Fig. |2t contribute to 
the disorder term, Fig. Due to the symmetry of the 
term in Eq. Q , the four possible types of diagrams 
in Fig. 12; contribute the same value. The result is 

x f dhdh j (2 7 r)' i (5(ki + k 2 + k 3 + k 4 ) 

x / dt 3 dt A J (27r) rf <5(k 5 + k 6 + k 7 + k 8 ) 

J ^kskekrks 

xlki + kal-'fe^Sj (ki,4i)Sj^(k 2 ,4i) 
x|k 5 + k 6 | y k 6s k Sv b r (k 5 ,t 3 )6i(k 8 ,4 4 ) 
xe-^-^0(t 3 - t 2 )n lt (k 6 )(2ir) d 5(k 6 + k 3 ) 
xe -^4(*2-t4) 0(t2 _ i 4 )77 lm (k4)(2^) rf <5(k 4 + k 7 ) 



nimn t-jji 



-K d 
2v 2 



—,S,ns5ltS„ 

a 



For the remaining contributions, we list the symmetry 
factor, terms in braces in the integrand of Eq. I)42|) that 
change, and the final contribution in brackets that change. 
The contributions are shown in Table where again the 
dependence on time has been suppressed. The terms e~ uki (* ! 
and e ~ vk i( t3 ~ t2 \ which have the same momentum argu- 
ments as the two II (kj) terms, have been suppressed. The 



d(d+ 2 
2 

'd(d+ 2) 
2 



d(d+2)(d + 'i) 

f A 

x / dqq d - y - 3 

J A/a 



-Ml 



^ijk u rst 



-*3 



X 



/d4id4 2 f (2 7 r)^(k 1 +k2+k3+k 4 ) 

J Jkik2k 3 k4 



|ki + k 2 |- y k 2j k 4u 



xSi (ki, ti)S£(k 2 , h)b r (k 3 , t 2 )bi(k 4 , t 2 ) (43) 
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It is clear that to evaluate this expression, we need to 
evaluate a term such as 

fci&.Mkl, *l)6fc (k2, h)br0<3, *2)Si(k4, *2) 

(44) 



XD ijk D rst [OltOmsOn 



Fourteen other terms need to be evaluated in order to cal- 
culate the total contribution from Eq. I|4H|I. These terms 
each contribute in a form that can be cast as a contribu- 
tion to D a and Dp in Eq. J5J. Shown in Tableware the 
terms and their contributions. Summing all these contri- 
butions to D a and Dp, we get the third and fourth flow 
equations of Eq. . 
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Table 1. The matrix M in the calculation of the intermittency exponent £4 via Eg, 1381 



M 



4(d 3 + 5d 2 + 2d - 6) 8(d + 3) 2 

2(d + 2) 2 2(d 3 + 5d 2 -12) 

4(d 3 + 5d 2 + 2d - 4) 8(d 2 +6d+ 10) 



2d(d + 2) 
2d(d + 2) 
8 

-2(d + 2) 



2(d 3 + 5d 2 - 2d - 16) 
2(d 3 + 5d 2 - 2d - 16) 



-4(d + 2) 
2(d 3 + 5d 2 + 2d) 
2d(d + 2) 
2d(d+ 2) 



16 

-8(d + 2) 
8(d 2 +6d+ 12) 



4(d 2 + 6d + 10) 



16 



d 3 +4d 2 -8d-24 2(d + 2) 2 4(d + 2) 2 
4(d 2 + 6d + 12) 16 32 



2(d 3 + 5d 2 - 6d - 24) d 3 + 4d 2 - 12d - 32 2d(d + 2) 4d(d + 2) 
2(d 3 + 5d 2 - 6d - 24) d 3 + 4d 2 - 12d - 32 2d(d + 2) 4d(d + 2) 



2(d 3 + 5d 2 + 2d - 4) 8(d 2 + 6d + 10) 

j3 



16 



-2(d + 2) (d + 2) 2 d 3 + 5d 2 -2d-16 

Table 2. Terms in the contributions from Fig. [5^,, via Eg. 1401 



-2(d + 2) 



-2(d + 2) -4(d + 2) , 



Integral 


Symmetry 


Integrand 


Result 


h 


2 


1^2^4,6, (ki)S+ 


(k 4 )i7 13 (k 2 )| 


^ a d(d+2) 




h 


2 


{k 2l k Al b i (ki)Sj- 


(k 4 ) J fT JJ (k 2 )} 




ra 
J 3 


1 


s kiik^jbk (ki)&i 


(k4)JT J fc(k 2 )| 


Do » 




2 d(d+2) 




T b 


1 


< k 2i k 4:) b k (k.3)bj 


(k 2 )77 lfc (k 4 )} 






2 d(d+2) 




J-i 


1 


{k 2% k A %j (ki)Sfc 


(k 4 )n Ifc (k 2 )| 





J b 

J-4 


1 


< k k 2i k4 j b i (k 3 )6 fe 


(k 2 )ZZ,- fe (k 4 )} 





ra 


1 x 2 


(k 2i k 4j b J (ki)S+ 


(k 4 )i7 fcfe (k 2 ) I 







1 x 2 


(k 2i k 4j b k (k 3 )Sfe 


(k 2 )i7y(k4)} 





T a 

6 


1 


< k 2i kAjbi (ki)6 fc 


(k 4 )i7, fc (k 2 )| 


2 d(d+ 2 ) 




T b 

J-a 


1 


•i k2 i k4 j b j (k 3 )6 fc 


(k 2 )77 lfc (k 4 )} 


Dp y 




2 d(d+2) 




h 


2 


< k 2l k4 t bi (ki)bj 


(k 4 )i7 l3 (k 2 )| 


d(d+2) 




, a 
J « 


1 x 2 


(k 2i k 4j bi (ki)6fc 


(k 4 )77, fc (k 2 )| 


d(d+ 2 ) 




T b 

■<8 


1x2 


f - x 

< k 2i k 4j b k (kzPj 


(k 2 )i7 lfc (k 4 )| 


U P d(d+ 2 ) 
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Table 3. Terms in the contributions from Fig. 12b, via Ea. 1421 



Integral Symmetry Integrand Result 



Jl 


2 


{ k 2k k 4k bf (kx )bf (k 4 )6^ (k 7 )iT j; (k 2 )J7 im (k 6 ) } 


p (d+l) 
^ Q d(d+2) 


J 2 


2 


{ k 2k k 4k bf (kx)^ (k 4 )fe^ (k 7 )i7, ; (k 2 )i7 mj (k 6 ) | 


[-^^; 


Js 


1 


{ k 2l k 4j if (kj )hi (k 4 )6^ (k T ) n jX (k 2 )n km (k 6 ) | 


+ Cc 1 

~ 2 d(d+2) 


T h 
■':', 


1 


|fe 2i fc 4 ife (k 3 )6 J ± (k 2 )^(k7)n«(k4)iT fcm (k 6 )| 


" + D a 1 
' 2 d(d+2) 


T" 
J 4 


1 


|fc2 I fc4i J ± (ki)^(k 4 )^(k 7 )i7 ifc (k 2 )iT lm (k 6 )| 





■h 


1 


{ fe 2j bf (k 3 )6^ (k 2 )6^ (k 7 )77 fc; (k 4 )i7 mj (k 6 ) I 





Js 


1x2 


|fc 2l fc 43 S J ± (ki)6 l ± (k 4 )^(k 7 )i7 fci (k 2 )i7 mfc (k6)| 





Js 


1x2 


[ feu fc 4j 6fc (k 3 )6^ (k 2 )6^ (k 7 )iT« (k 4 )77 mj (k 6 ) | 





■hi 


1 


{ k 2% k 4j bf (kr)S^ (k 4 )6^ (k 7 )77 fc; (k 2 )77 m3 (k 6 ) | 


" | D H 1 

"r" 2 d(d+2) 


T h 
■hi 


1 


|fc 2s fc 4 i J ± (k 3 )^(k 2 )^(k 7 )i7 fe; (k 4 )iT lm (k 6 )| 


[ , D H 1 
~ 2 d(d+2) 


Jr 


2 


{ k 4k bf (ki )S+ (k 4 )6^ (k 7 )iT« (k 2 )i7 jm (kfi) } 


f Da 

d(d+2) 


Js 


1x2 


| fc 2j fe^ Si (ki)S^ (k 4 )&^ (k 7 )^ii (k 2 )i7 fcm (k e ) | 


[ +I) /3d(d+2) 


Js 


1x2 


{fca^S* (k 3 )6 J ± (k 2 )^(k 7 )i7 fc; (k 4 )i7 ml (k 6 )| 


+I) /3 d(d+2) 
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Table 4. Terms in the contributions from Fig.|2t, via Ea.1441 



Term 


Contribution to D a 


WUIILI lULlLrlUll LL> i_y ^3 




A, A A 

u lt u ms u nu 


wl 


+ 4D a £> /3 


4 7~) 2 .L ('/i -L ^ n n,, _i_ ^ n 2 

^i^ Q -h 1^4 -|- a)U a Up -\- oijq 




A, A A 

u lt u mn u su 







(a _i_ n 2 _l ffi _l 0/?^ n n , j_ on 2 
^4 -+- za)U a -f- t^D -+- za)U a iy/3 + ^^73 




A A A 

u lt u mu u ns 


2Dl 


+ 4D a D /3 


4.L/ a -f- (4 -f- a)U a Ufj + O 




AAA 







/I 7~1 2 _i_ /'A _i_ OW^ 7~) 7~) _i_ fO J_ 7~) 2 
4i^ Q + (0 + Za)U a Ug + 4- ZajJJp 




A, A A 

u lm, u nt u su 







^4 -+- za)U a ■+■ (^4 -f- oa -+- a )U a u @ -+- yi - 


r a) -Up 


AAA 

Vlm.utu'-'ns 







(4 + 2d)D" + (4 + 3d + d 2 )D a D p + (1 - 


Vd)D 2 


A, A A 


21)2 




4D^ + (4 + d)D a D fj + 3D 2 3 




A A A 
tJ ls u mn vtu 







(4 + 2d)££ + (6 + 2d)D a D p + 2D} 




A A A 




+ 4£> a £)^ 


Wi + (4 + d)D a D + ZD% 











(4 + 2d)Z^ + (4 + 3d + d 2 )D a Dp + (1 - 


Vd)D 2 









(4 + 2d)L£ + (4 + 3d + d 2 )D a D p + (1 + d)L^ 


Sln$muSst 







4L£ + (6 + 2d)D a D p + (2 + 2d)Z^ 






2D 2 a 


+ 4L> a D /3 


4l£ + ( 4 + 2d)D a D + 2D} 




&lu $mn $st 







4D 2 a + 8D a D + AD% 








+ 4L> a -D/3 


Wl + (4 + 2d)D a Dp+2Dl 





